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ABSTRACT 

Using evolutionary population synthesis we present integrated colours, integrated spec- 
tral energy distributions and absorption-line indices denned by the Lick Observatory 
image dissector scanner (referred to as Lick/IDS) system, for an extensive set of in- 
stantaneous burst binary stellar populations with and without binary interactions. 
The ages of the populations are in the range 1 — 15 Gyr and the metallicities are in 
the range 0.0001 — 0.03. By comparing the results for populations with and without 
binary interactions we show that the inclusion of binary interactions makes the inte- 
grated U — B, B — V, V — R and R — I colours and all Lick/IDS spectral absorption 
indices (except for H^) substantially smaller. In other words binary evolution makes 
a population appear bluer. This effect raises the derived age and metallicity of the 
population. 

We calculate several sets of additional solar-metallicity binary stellar populations 
to explore the influence of input parameters to the binary evolution algorithm (the 
common-envelope ejection efficiency and the stellar wind mass-loss rate) on the result- 
ing integrated colours. We also look at the dependence on the choice of distribution 
functions used to generate the initial binary population. The results show that varia- 
tions in the choice of input model parameters and distributions can significantly affect 
the results. However, comparing the discrepancies that exist between the colours of 
various models, we find that the differences are less than those produced between the 
models with and those without binary interactions. Therefore it is very necessary to 
consider binary interactions in order to draw accurate conclusions from evolutionary 
population synthesis work. 
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1 INTRODUCTION 

Amongst the distinct methods available for the study of the 
integrated light of stellar populations it is the evolution- 
ary population synthesis (EPS) method (first introduced 
by Tinsley 1968) that offers the most direct approach for 
modelling galaxies. Remarkable progress in this field has 
been made during the past decades but the majority of EPS 



gle stars (Bressan, Chiosi & Fagotto 


1 1994 IWorthevlll994l 


Vazdekis ct al. 


1996: iKurth. Alvensleben & 


Frickell 199S(h 



Both observation and theory tell us that binary stars 
play a very important role in the evolution of stellar popu- 
lation. First, observations show that upwards of 50% of the 
stars populating galaxies are e xpected to be in binary or 
higher-order multip le systems .Duauennov fe Mavorlll991_ 
iRichichi et al] |l994. for example). Secondly, binary evolu- 
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tion, if the component stars are close enough to exchange 
mass, can drastically alter the evolution path of a star as 
expected from single star evolution. Moreover, binary inter- 
actions can also create some important classes of o bjects 
such as blue stragglers fBSs: |p"ols fc Marinusl |l994}h and 
subdwarf B stars (sd Bs, also referred as ext reme horizontal 
branch [EHB] stars: lHan et alJl2002. 1200.'-) . These objects 
produced by binary evolution channels, can significantly af- 
fect the integrated spectral energy distribution (ISED) of 
a population from ultraviolet (UV) to radio ranges. There- 
fore it is necessary to include binary stars in EPS mod- 
els. However, in the few EPS studies to date that have ac- 
counted for binary evolution, some only included specialised 
classes of binaries (for example, massive clo se binary evo- 
lution: IVan Bever fc Vanbeverenlll99Sl l200ff) . while others 
were limited to studyi n g specialised stellar populations (e.g . 

Cervino. Mas-Hesse fc Kuntb] 119971 : 
The EPS study of 



Pols & Marinus 1994: 



Van Bever fc Vanheverenl ll 99TE 

Zhang et alJ <l2004bt hereafter Paper I) took into account 
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various known classes of binary stars and investigated gen- 
eral populations but they only considered solar metallicity 
populations. This paper expands that study by including a 
metallicity dependence and by investigating the effects of 
binary interactions and model input parameters on the re- 
sults. Furthermore, by using 1 x 10 6 instead of 2 x 10 s binary 
systems the error in the calculations is reduced and the solar 
metallicity results of Paper I are superceded. 

In this paper we assume that all stars are born in bi- 
naries and born at the same time, i.e. an instantaneous bi- 
nary stellar population (BSP). EPS models of instantaneous 
burst BSPs require four key ingredients: (i) a library of evo- 
lutionary tracks (including binary stars and single stars) 
used to calculate isochrones in the colour-magnitude dia- 
gram (CMD) ; (ii) a library of stellar spectra adopted in order 
to derive the ISED, or magnitudes and colours, in suitable 
passbands; (iii) a method to transform spectral information 
to absorption-line strengths, for example, the ap proach of 
empirical fitting polynomials llWorthev et al Il994l) ; and (iv) 
assumptions regarding the various distributions required for 
initialization of the binary population, such as the initial 
mass function (IMF) of the primaries and the distribution 
of orbital separations. 

The outline of the paper is as follows: we describe our 
EPS models and algorithm in Section 2; our results are pre- 
sented in Section 3 and in Section 4 we investigate the in- 
fluences of binary interactions, model input parameters and 
distribution functions on the results; and then finally, in Sec- 
tion 5, we give our conclusions. 



2 MODEL DESCRIPTION 

The EP S code used here was develop ed by Zhang and col- 
leagues iZhang et al] I2002I . l2004al lbr) and is based on the 
most updated input physics. Here we briefly describe the 
important features and components of the model as well as 
providing details of the algorithm for computing a BSP. 



2.1 Initialization of the binary population 

To investigate the effect of binary interactions on the EPS 
study of stellar populations we first need to construct the in- 
stantaneous burst BSPs. The evolution path taken by a star 
in a binary depends critically on the mass of the compan- 
ion star and the orbital parameters. As such, binary evo- 
lution can be very complex and it is difficult to estimate 
the evolutionary timescales of the two stars without actu- 
ally evolving the system. Therefore, the traditional method 
of constructing populations used in the EPS study of single 
stellar populations, where only the IMF need be considered 
and the evolutionary timescale of a star is set as soon as 
its mass is set, confronts difficulty in the study of BSPs. 
Here we need to use a Monte Carlo process, which utilises a 
random-number generator in combination with distribution 
functions, to generate the initial conditions of a set of bi- 
naries which then need to be evolved using an appropriate 
binary evolution model. 

To generate a BSP the following input distributions are 
required to define the initial state of each binary: (i) the 
IMF of the primaries, which gives the relative number of 
the primaries in the mass range M — * M + AM; (ii) the 



secondary-mass distribution; (iii) the distribution of orbital 
separations (or periods); and (iv) the eccentricity distribu- 
tion. For distributions for which the majority of studies have 
reached fair agreement, such as the primary-mass and or- 
bital separation distributions, we use a certain form while 
for those that are less constrained (secondary-mass and ec- 
centricity) we use several reasonable assumptions. We also 
need to set the lower and upper mass cut-offs M\ and M u 
to the mass distributions and assign a metallicity Z to the 
stars. 

For each binary system the initial mass o f the primary is 
chosen from the ap proximation to the IMF oflMill er fc Scale] 
lll979f) as given bv lEggleton. Fitchett fc Toutl <1989t here- 
inafter EFT), 



Mi 



Q.19X 



(1 -X) ' 75 + 0.032(1 - X) c 



(1) 



where X is a random variable uniformly distributed in the 
range [0,1], and Mi is the primary mass in units of Mq. 

In choosing the initial mass of the secondary star we 
can assume that the masses of the component stars are ei- 
ther correlated or that they are independent. For the corre- 
lated case, the secondary-mass distribution depends on the 
primary-mass (as set by equation and the mass-ratio, 
q, distribution. The form of the latter is somewhat uncer- 
tain and a matter for debate. In this study we consider two 
versions of the o distribution, one is a uniform distribut ion 
(EFT 1989: lMazeh et al.lll992l : iGoldberg fc Mazehlll994D . 



n(q) 



< q < 1, 



(2) 



where q = M2/M1, and the other is a thermal distribution 
<Han. Podsiadlowski fc Eggletodl7995l) . 



n(q) = 2q, < q < 1. 



(•3) 



In the uncorrelated case, the secondary mass is chosen in- 
dependently from the same IMF as the primary (equation 

The distribution of orbital separations is taken as con- 
stant in log a (where a is the separation) for wide binaries 
and falls off smoothly at close separations: 



i(a) = { ascpl 



(a/a ) 



a < ao, 

ao < a < en, 



(4) 



where a scp w 0.070, ao 



10R Q , ai = 5.75 x 10 6 R Q and 
m ~ 1.2. This distribution implies that there are equal num- 
bers of wide binary systems per logarithmic interval, and 
that approximately 50% of the binary systems have orbital 
periods less than 100 yr. The value of 100 yr can be viewed 
as the upper limit for interaction between the component 
stars - if the period is longer than 100 yr the evolution of 
the stars does not differ from that of two independent sin- 
gle stars. This fraction of 50% for binaries with a period 
less than 100 yr is a typical value for the Galaxy, resulting 
in ~ 10% of the binaries experiencing Roche lobe overflow 
(RLOF) during the past 13 Gyr. 

In order to investigate the effect of eccentricity we con- 
sider two eccentricity distributions: (i) all binaries are ini- 
tially circular, i.e. 



e = 0, 



(5) 



and (ii) all binaries are formed in eccentric orbits where the 
initial eccentricity distribution satisfies a uniform form, i.e., 
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x. 



(6) 



Under this second assumption large eccentricities in short- 
period orbits are to be excluded on the basis that the stars 
will crash into each other at periastron. 

Equations - © can be used to set the initial state 
of a binary system: the masses of the component stars, Mi 
and M2, separation, a, and eccentricity, e, of the orbit. 



2.2 Binary evolution model and input parameters 

To describe the evolution of a binary we use the rapid bi- 
nary s tar evolution (BSE) algorithm of Hu rley. Tout fc Polsl 
(2002). The BSE algorithm provides the stellar luminosity 
L, effective temperature T e g, radius R, current mass M and 
the ratio of radius to Roche-lobe radius R/Rl for the com- 
ponent stars, as well as the period P, separation a and ec- 
centricity e for a binary system. It is valid for component 
star masses in the range 0.1 < Mi, M2 < lOOM©, metallicity 
0.0001 < Z < 0.03, and eccentricity 0.0 < e < 1.0. The BSE 
algorithm includes the single star evol ution (SSE) pa ckage 
of ana lytic formulae as presented by 'Hurley. Pols fc Toutl 
(2000) in its entirety. In fact, for orbits that are wide enough 
that mass exchange between the component stars does not 
take place, the evolutionary parameters of the stars are iden- 
tical to that given by the SSE package. The SSE package 
comprises a set of analytic evolut ion functions fitted to the 
model tracks of iPols et alJ Jl998t) . Detailed descri ptions of 
the stellar evolutionary models of lPols et alJ (I1998T) and the 
SSE package have been presented previously in Zhang et 
al. (2002: see Sec. 2), so we do not discuss them here. In 
addition to all aspects of single star evolution, the BSE al- 
gorithm models processes such as mass transfer, mass ac- 
cretion, common- envelope (CE) evolution, collisions, super- 
nova kicks, tidal evolution, and all angular momentum loss 
mechanisms. This is done mostly by using a prescription (or 
recipe) based approach. 

For the BSE code there are several important input 
parameters that require mention: 

• (i) the efficiency of CE ejection qce in the CE evolution 
model denotes the fraction of the orbital energy that is trans- 
ferred to the envelope and is use d to overcome the b inding 
energy. This efficiency is defined jlben fc Livid ll993) by 



A-Ebin 

AE OI ] 



(7) 



where A_E or b is the change in the orbital energy of the binary 
between the initial and final state of the spiraling-in process 
and AiSbind is the energy added to the binding energy of the 
envelope. If the envelope is ejected before complete spiral- 
in then a close binary is the result, otherwise the two stars 
will coalesce. This parameter is crucial in understanding the 
evolution of populations of binary systems because it deter- 
mines the outcome of a CE interaction. In our BSP study 
we vary it over a reasonable range (1.0 — 3.0) to investigate 
its effects. 

• (ii) the c oefficient 77 for Reimers' wind mass-loss 
jReimerslll975h . 



M R = 4x 10~ 13 r, — M Q yr-\ 



(8) 



solar units. Equation JSJ gives the wind mass-loss rate for 
intermediate and low-mass stars on the giant branch and 
beyond. In this study we vary r\ from 0.5 to 0.0. 

• (iii) the tidal enhancement parameter B, which ap- 
pe ars in the formula of t idally enhanced mass loss given 
bv lTout fc EggletorJ dl988l> . 



M = Mr 



1 + B w max(i,-^-) 



(9) 



where Mr is the Reimers' rate (see equation and 7?l 
is the Roche lobe radius. Equation © is used to increase 
the incidence of RS CVn binaries. In this study we do not 
consider tidally enhanced mass loss, i.e., B — 0.0. 

2.3 Stellar spectra and absorption line indices 

The BaSeL-2 .0 stellar spectra library of iLeieune et alJ 
1 1991 r 1998T) provides an extensive and homogeneous grid 
of low-resolution theoretical flux distributions in the range 
of 9. 1 - 160000 nm, and synthetic UBVRIJHKLM colours for 
a large range of stellar parameters: 2000 < T c ff/K < 50000, 
-1.02 < log 3 < 5.50, and +1.0 < [Fe/H] < -5.0 (where 
g denotes surface gravity). For this library correction func- 
tions have been calculated for each value of the T e g and for 
each wavelength in order to yield synthetic UBVRIJHKLM 
colours matching the empirical colour- T c g calibrations de- 
rived from observations at solar metallicity. Semi-empirical 
calibrations for non-solar abundances ([Fe/H] = —3.5 to 
+ 1.0) have also been established for this version of the li- 
brary. After correction the most important systematic differ- 
ences existing between the original model spectra and the 
observations are eliminated. Furthermore, synthetic UBV 
and Washington ultraviolet excesses S 



(U-B)> O(C-M) 



and 



where stellar luminosity, L, radius, R, and mass M, are all in 



Src-Tx), obtained from the original model spectra of giants 
and dwarfs, are in excellent agreement with the empirical 
metal-abundance calibrations. 

The empirical fitting functions of lWorthev et all (11994) 
give absorption-line indices defined by the Lick Observa- 
tory image dissector scanner (Lick/IDS) system as a func- 
tion of T c g, logg, and metallicity [Fe/H]. The effective tem- 
perature spans a range of 2100 < T e a/K < 11000 and the 
metallicity is in the range —1.0 < [Fe/H] < +0.5. The in- 
dices in the Lick system were extracted from the spectra of 
460 stars obtained between 1972 and 1984 using the red- 
sensitive IDS and Cassegrain spectrograph on the 3m Shane 
telescope at Lick Observatory. The spectra cover the range 
4000 - 6400 A, with a resolution of ~ 8 A JWorthev et al] 
11994) . A more detailed description of the Lick/IDS spectral 
absorption-feature indices has been presented in Zhang et 
al. (2004a: see Sec. 2). 



2.4 The EPS algorithm 

After choosing a seed for the random number generator we 
use the Monte Carlo eqns Q - © to produce a popula- 
tion of 1 x 10 6 binary systems. In this work the lower and 
upper mass cut-offs Mi and M u are taken as 0.1 Mq and 
100 q, respectively. The relative age, r, of the BSP is as- 
signed within the range of 1 — 15 Gyr and the metallicity is 
chosen within the limits 0.0001 < Z < 0.03. We can then 
use the BSE algorithm to evolve each binary in the BSP 
to an age of r which gives us evolutionary parameters such 
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as L, T G ff, 7? and M for the component stars. Next we use 
the BaSeL-2.0 stellar spectra to transform the evolutionary 
parameters to colours and stellar flux and t hen u se the em- 
pirical fitting functions of IWorthev et aD (|j994) to derive 
spectral absorption feature indices in the Lick/IDS system. 
Finally, by the following equations HUH - 1131 we can obtain 
the integrated colours, monochromatic flux and absorption 
feature indices for an instantaneous BSP of a particular age 
and metallicity. We note that the common-envelope phase 
of binary evolution is assumed to be instantaneous in the 
BSE algorithm and thus it is not a factor when calculating 
the flux. Other aspects of binary evolution that may affect 
the observed flux, such as the presence of an accretion disc 
or novae eruptions, are also not taken into account, i.e. we 
only consider the stellar parameters. 

In the following equations, a parameter identified by a 
capital letter on the left-hand side represents the integrated 
BSP, while the corresponding parameter in minuscule on the 
right-hand side is for the k— th star. The integrated colour 
is expressed by 



(d - Cj) T ,z 



-2.5 log 



(10) 



where a and Cj are the i-th and j-th magnitude of the fc-th 
star. 

The integrated monochromatic flux of a BSP is defined 

as 



Fx- 



(ii) 



where fx is the SED of the fc-th star. 

The integrated absorption feature index of the Lick/IDS 
system is a flux-weighted one. For the i— th atomic absorp- 
tion line, it is expressed in equivalent width (W, in A), 



(12) 



where w\ is the equivalent width of the i— th index of the 
fc-th star, and fi t cx is the continuum flux at the midpoint 
of the i— th 'feature' passband; and for the i— th molecular 
line, the feature index is expressed in magnitude, 



Ci. 



-2.5 log 



ELi io- 



(13) 



where Cj is the magnitude of the i— th index of the fc-th star 
(as in equation 1101 . 
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Figure 1. The effects of binary interactions on U — B colour for 
BSPs of various metallicity. The symbols linked by line are for 
BSPs with binary interactions (Model A') and those without a 
line are without binary interactions (Model B'). Different symbols 
denote different metallicities, from top to bottom, the metallicity 
Z is 0.03, 0.02, 0.01, 0.004, 0.001, 0.0003 and 0.0001, respectively. 



I 




Age (Gyr) 

Figure 2. Similar to Fig.Qbut for B — V colour. 



3 RESULTS 

In this part we present integrated colours, ISEDs at interme- 
diate resolution (10 A in the UV and 20 A in the visible) 
and Lick/IDS absorption feature indices for instantaneous 
burst BSPs with and without binary interactions over a 
large range of age and metallicity: 1 < r < 15 Gyr and 
-2.3 < [Fc/H] < +0.2. For each model, a total of 1.0 x 10 6 
binaries are evolved according to the algorithm given in the 
previous section. 



3.1 Colours 

In Table □ we present U - B, B - V, V - R and V - I 
colours for the set of BSPs that we call Model A'. In this 
model binary interactions are included by utilizing the BSE 
algorithm, the initial mass of the secondary M2 is assumed 
to be correlated with Mi and is obtained from a uniform 
mass-ratio distribution (see equation eccentric orbits are 
allowed for the binary systems and a uniform eccentricity 
distribution is adopted (see equation |SJ , the CE ejection 
efficiency acE is set to 1.0, and the Reimers' wind mass-loss 
coefficient rj is 0.3. The model is repeated for seven distinct 
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Table 1. The integrated colours for Model A'. 



Age 
(Gyr) 


1.00 


2.00 


3.00 


4.00 


5.00 


6.00 


7.00 


8.00 


9.00 


10.00 


11.00 


12.00 


13.00 


14.00 


15.00 


















Z = 


0.0001 
















U - 


B 


-0.021 


-0.001 


0.009 


0.018 


-0.009 


-0.006 


-0.007 


0.005 


0.007 


-0.012 


0.007 


-0.004 


0.002 


-0.010 


-0.018 


B - 


V 


0.266 


0.404 


0.484 


0.547 


0.515 


0.554 


0.553 


0.569 


0.579 


0.534 


0.573 


0.574 


0.580 


0.587 


0.572 


V - 


R 


0.274 


0.339 


0.382 


0.408 


0.366 


0.386 


0.389 


0.401 


0.404 


0.373 


0.400 


0.403 


0.410 


0.411 


0.399 


V - 


- 1 


0.647 


0.755 


0.829 


0.870 


0.781 


0.816 


0.824 


0.846 


0.853 


0.790 


0.841 


0.849 


0.861 


0.861 


0.838 


















Z = 


0003 
















u - 


B 


-0.063 


0.011 


0.015 


0.025 


0.011 


-0.008 


-0.011 


-0.030 


-0.010 


0.004 


0.007 


0.010 


0.046 


0.030 


0.032 


B - 


V 


0.295 


0.431 


0.518 


0.583 


0.593 


0.625 


0.631 


0.607 


0.654 


0.667 


0.666 


0.654 


0.691 


0.651 


0.646 


V - 


R 


0.271 


0.350 


0.395 


0.427 


0.407 


0.425 


0.431 


0.406 


0.435 


0.448 


0.441 


0.439 


0.470 


0.442 


0.437 


V - 


- 1 


0.627 


0.781 


0.855 


0.911 


0.856 


0.885 


0.899 


0.841 


0.898 


0.922 


0.909 


0.911 


0.970 


0.916 


0.904 


















Z = 


0.001 
















u - 


B 


-0.038 


0.045 


0.076 


0.065 


0.070 


0.068 


0.051 


0.050 


0.033 


0.045 


0.073 


0.038 


0.058 


0.062 


0.108 


B - 


V 


0.343 


0.463 


0.577 


0.632 


0.663 


0.679 


0.692 


0.712 


0.688 


0.692 


0.727 


0.677 


0.658 


0.678 


0.664 


V - 


R 


0.303 


0.360 


0.401 


0.437 


0.440 


0.445 


0.449 


0.470 


0.448 


0.448 


0.473 


0.440 


0.431 


0.441 


0.446 


V - 


- 1 


0.696 


0.821 


0.872 


0.952 


0.934 


0.930 


0.931 


0.974 


0.925 


0.920 


0.969 


0.908 


0.887 


0.908 


0.929 


















Z = 


0.004 
















u - 


B 


0.077 


0.116 


0.102 


0.134 


0.164 


0.198 


0.211 


0.192 


0.221 


0.242 


0.257 


0.253 


0.275 


0.275 


0.273 


B - 


V 


0.366 


0.544 


0.624 


0.692 


0.730 


0.775 


0.790 


0.772 


0.803 


0.824 


0.833 


0.819 


0.833 


0.844 


0.838 


V - 


R 


0.281 


0.356 


0.401 


0.430 


0.449 


0.471 


0.483 


0.476 


0.482 


0.495 


0.501 


0.507 


0.511 


0.514 


0.509 


V - 


- I 


0.689 


0.808 


0.891 


0.920 


0.959 


0.990 


1.034 


1.020 


1.001 


1.026 


1.049 


1.081 


1.065 


1.074 


1.055 


















Z = 


- 0.01 
















u - 


B 


0.139 


0.147 


0.198 


0.237 


0.275 


0.308 


0.305 


0.348 


0.372 


0.384 


0.394 


0.408 


0.433 


0.443 


0.458 


B - 


V 


n A a i 


U.D4Z 


fl 70Q 
U. I Zo 


n TVLT 
U. 1 o I 


n 7Q/1 
U. I o4 


n coo 
u.ozy 


U.ozD 


u.ooy 


U.ool 


U.ooD 


U.ooo 


U.oyo 


u.yuo 


n ni i 

u.yii 


n oi7 
u.y 1 1 


V- 


R 


0.295 


0.399 


0.439 


0.457 


0.471 


0.493 


0.491 


0.507 


0.518 


0.522 


0.521 


0.525 


0.529 


0.529 


0.537 


V- 


- 1 


0.704 


0.890 


0.936 


0.980 


0.991 


1.033 


1.024 


1.047 


1.067 


1.082 


1.080 


1.081 


1.091 


1.089 


1.111 


















Z = 


= 0.02 
















u- 


B 


0.167 


0.204 


0.289 


0.351 


0.396 


0.438 


0.487 


0.480 


0.512 


0.587 


0.587 


0.618 


0.639 


0.642 


0.705 


B - 


V 


0.525 


0.696 


0.788 


0.839 


0.878 


0.905 


0.938 


0.919 


0.960 


0.990 


0.998 


0.981 


1.022 


1.029 


1.037 


V- 


R 


0.322 


0.422 


0.469 


0.491 


0.510 


0.523 


0.541 


0.536 


0.554 


0.566 


0.570 


0.559 


0.577 


0.565 


0.586 


V- 


- 1 


0.704 


0.889 


0.985 


1.006 


1.052 


1.071 


1.109 


1.092 


1.135 


1.165 


1.168 


1.129 


1.166 


1.160 


1.186 


















Z = 


= 0.03 
















u- 


B 


0.183 


0.278 


0.355 


0.415 


0.470 


0.495 


0.518 


0.557 


0.648 
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values of Z: 0.03, 0.02, 0.01, 0.004, 0.001, 0.0003 and 0.0001, 
respectively. 

To investigate the effect of binary interactions on the re- 
sults, we also construct Model B', which differs from Model 
A' by neglecting all binary interactions, i.e. the component 
stars are evolved as if in isolation according to the SSE algo- 
rithm. In Figs0and|5|we give a comparison of the integrated 
U — B and B— V colours for Models A' (generated from Table 
and B'. We see that the inclusion of binary interactions 
makes the integrated U — B and B — V colours of various 
metallicity BSPs bluer for all instances. 

In Figs |3] and |3] we give the fractional contribution of 
different evolutionary stages to the total flux for solar metal- 
licity 1— and 13— Gyr BSPs, respectively. The top panels of 
Figs |3] and 0] are for Model A', while the bottom panels are 
for Model B'. In each of them various abbreviations are used 
to denote the evolution phases. They are as follows: 'MS' 
stands for main-sequence stars, which are divided into two 
phases to distinguish deeply or fully convective low-mass 
stars (M<0.7Mq, therefore Low-MS) and stars of higher 
mass with little or no convective envelope (M>O.7M0); 
'HG' stands for Hertzsprung gap; 'GB' stands for the first gi- 
ant branch; 'CHeB' stands for core helium burning; 'EAGB' 



stands for early asymptotic giant branch stars; 'PEAGB' 
(post EAGB) refers to those phases beyond the 'EAGB' - in- 
cluding thermally pulsing giant branch/proto planetary neb- 
ula/planetary nebula (TPAGB/PPN/PN); 'HeMS' stands 
for helium star main-sequence; 'HeHG' stands for helium 
star Hertzsprung gap; 'HeGB' stands for helium star giant 
branch; and 'HeWD' and 'CO WD' stand for helium and car- 
bon/oxygen white dwarfs, respectively. Correspondingly, in 
Figs 03 and |5] we plot theoretical isochrones of solar metal- 
licity 1— and 13— Gyr BSPs. Note that all stars except for 
Low-MS stars are included in the isochrones. 

Comparing the contribution of different evolutionary 
stages to the total flux between Model A' and Model B' 
for solar metallicity 1- Gyr BSPs (Fig. we find that the 
flux in the UBVRI passbands is dominated by the same 
evolutionary stages regardless of model. This is also true for 
the 13 Gyr solar metallicity BSPs (Fig. 2J. For the young 
BSPs the UBVRI flux is dominated by MS stars with mass 
M>O.7M , CHeB stars and cooler PEAGB (in fact the 
TPAGB) stars. MS and CHeB stars also dominate for the 
old BSPs but here the GB stars make a major contribu- 
tion as well. Therefore, the differences in the integrated 
U-B,B-V,V-R and V - 1 colours between Model A' 
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Figure 3. The fractional contributions of different evolutionary 
stages to the total flux are shown for solar metallicity BSPs at an 
age of 1 Gyr. Top and bottom panels are for Models A' and B', 
respectively. The abbreviations are explained in the text. 



and Model B' will originate from differences in the distri- 
bution of these classes of stars which dominate the total 
UBVRI flux. From Figs|S]and|H]we see clearly that the dis- 
tribution of stars in the Hertzsprung-Russell (HR) diagram 
for Model A' is significantly different from that of Model B'. 
The distribution of the stars is more dispersed for Model A' 
in comparison with Model B' and, as you would expect, BS 
stars are produced only by Model A'. Moreover, for young 
BSPs, Fig. ©shows that Model A' produces a lot of scatter 
around the MS as well as hotter CHeB stars compared to 
Model B'. For old BSPs, Fig. ©shows that Model A' also 
produces MS, GB and CHeB stars which differ from the cor- 
responding class of stars in Model B'. The net effect of the 
difference in the distribution of these classes of stars is to 
make the populations hotter and therefore Model A' appears 
younger when compared with Model B'. 

So far we have neglected to discuss redder colours (such 
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Figure 4. Similar to Fig. but for an age of 13 Gyr. 



as, J — K, H — K, etc) and the reason is their relatively 
large errors. From the top and bottom panels of Fig. © we 
see that cooler PEAGB stars {T e g < 3.54, see Fig. © give 
the maximal contribution to the total flux in those redder 
passbands for young BSPs. The evolutionary timescale of 
these PEAGB stars is so short compared with that of MS 
stars that there are few PEAGB stars in our 1 x 10 6 binary 
system sample. Therefore, for young BSPs large error can 
result from fluctuations in the small number of PEAGB stars 
and thus, the evolutionary curves of these colours are not 
smooth. For old BSPs, from the top and bottom panels of 
Fig. 0| we see that GB stars contribute about 50% of the 
total flux in the redder wavelength range: in fact, it is GB 
stars with temperature T e g < 3.54 that give the maximal 
contribution. Although the evolutionary timescale of these 
stars is longer than that of young PEAGB stars it is still 
relatively short compared to the age of the population and 
once again a large error can be introduced in the red flux 
range. The reason we focus on the bluer colours is that the 
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Figure 5. Theoretical isochrones for solar-metallicity instanta- 
neous burst BSPs at an age of 1 Gyr. Open circles for Model B', 
dots for Models A' (black for the primary, blue for the secondary). 
For the sake of clarity only 1 X 10 5 binary systems are included 
and those MS stars with mass M <0.7Mq are removed. All bi- 
naries are assumed to be resolved. 
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Figure 6. Similar to Fig.^but for an age of 13 Gyr. 

total flux in the UVBRI passbands is dominated by stars 
which have long evolutionary timescales and the error here 
is comparitively small. 

3.2 The integrated spectral energy distribution 

In Figs □ and |H1 we give the variation of the intermediate- 
resolution ISED with age and metallicity over a wide wave- 
length range, 2.7 < log A/ A < 5.2, for Model A'. The flux 
is expressed in magnitudes and is normalized to zero at 2.2 
/jm. Figs|2land|S|show that the effects of age and metallicity 
on the ISEDs are similar, i.e. the ISEDs tend to be redder 
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Figure 7. The integrated spectral energy distributions as a func- 
tion of age for solar-metallicity BSPs with binary interactions 
(Model A'). From top to bottom, the ages are r = 1,2,4,7,10 
and 15 Gyr, respectively. The flux is expressed in magnitudes 
and is normalized to zero at 2.2 /im. 
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Figure 8. The integrated spectral energy distributions as a func- 
tion of metallicity for BSPs with binary interactions (Model A') at 
age t = 10 Gyr. From bottom to top, the metallicity Z is 0.0001, 
0.0003, 0.001, 0.004, 0.01, 0.02 and 0.03, respectively. The flux is 
also expressed in magnitudes and is normalized to zero at 2.2 (im. 

with increasing age and metallicity in the wavelength region 
3.3 < log A/A < 4.2. 

In Fig. El we compare the ISEDs of Model A' with 
those of Model B' for solar-metallicity BSPs at ages r = 
1,2,4,7,10 and 15 Gyr (note that the ISEDs for Model A' 
have also been given in Fig. 0. In Fig. [5] for Model A' ar- 
rows mark the blue end (Ai, T ) at which MS stars with mass 
M > O.7M give a 10% contribution to the total flux. For 
Model B' we do not give the fluxes shorter than the wave- 
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Figure 9. The integrated spectral energy distributions for so- 
lar metallicity instantaneous BSPs with (full line) and without 
(dashed line) binary interactions at ages r = 1, 2, 4, 7, 10 and 15 
Gyr (from top to bottom, respectively). For Models A' arrows 
mark the blue end at which MS stars with mass M < O.7M0 con- 
tribute 10% of the ISED, and for Model B' we do not plot flux 
shorter than this point. Note for young BSPs the blue end where 
MS stars with mass M<O.7M0 give a 10% contribution is sig- 
nificantly different for Models A' and B', while for old BSPs this 
discrepancy is insignificant. For the sake of clarity the fluxes at 
ages of 2, 4, 7, 10, 15 Gyr are shifted downwards by an mount of 
-4.0. 



length A2, T (defined as for Ai jT but for Model B'). Note that 
for young BSPs Ai jT of Model A' is significantly different 
from A 2 , T of Model B', with Ai, T < A 2 , r , while for old BSPs 
we find Ai, T ~ A2, T - The differences between Ai jT and \2, T 
are explained by the absence of BS stars in Model B', while 
in Model A' the number of hot BS stars at early ages is 
greater than at old ages. Fig. [5] shows significant disagree- 
ment of the fSED between Model A' and Model B' in the 
UV and far-UV regions (log A/A < 3.4): the ISED for Model 
A' is bluer than that for Model B' at all ages. Comparing 
this discrepancy for the young and old BSPs it seems that it 
decreases with age: at r = 1 Gyr the discrepancy reaches to 
~ 5.0 mag at log A/A = 3.1 while at r = 15 Gyr it decrease 
to ~ 2.2 mag at log A/A = 3.3. 

To aid investigation of the differences in the ISED in 
the UV and far-UV regions between Model A' and Model 
B' we have included arrows to mark the wavelength corre- 
sponding to Ai jT (left) and \2, T (right) in Figs |3] and 0] For 
young r = 1 Gyr BSPs, Fig.EJshows that the ISED in the re- 
gion of A2,t < log A/A < 3.4 is dominated by MS stars with 
mass M >O.7M for both Model A' and Model B', and the 
fractional contribution of these MS stars drops from 0.93 to 
0.1 rapidly at A2, T for Model B'. So the discrepancy of the 
ISED in this region is mainly introduced by the difference 
in the distribution of MS stars with mass M > O.7M0 (see 
Fig.®. In the region Ai, T < A < A2, T the total flux is mainly 
dominated by hotter PEAGB and/or CHeB stars for Model 
B', while still by MS stars with mass M > O.7M for Model 
A', so the discrepancy in this region is mainly caused by the 



Table 2. Model Parameters. 
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existence in Model A' of more hot BSs. This analysis of the 
t = 1 Gyr BSP is also suitable for the r = 2, 4, 7, 10 Gyr 
BSPs. For the old r — 15 Gyr BSP, Ai iT is almost equal to 
A2,t and in this case the discrepancy of the ISED in the 
UV and far-UV regions is caused mainly by the difference in 
the distributions of MS stars with mass M >O.7M , hotter 
PEAGB and/or CHeB stars (as exhibited in Fig.®. 

For the young BSPs (r = 1 and 2 Gyr) we also see 
some evidence in Fig. of a disagreement of the ISEDs 
for Models A' and B' in the visible and infrared regions 
(3.5 < log A/A < 4.0). It appears that Model A' exhibits a 
bluer continuum than Model B' in this range and we be- 
lieve this is owing to fluctuations introduced by the small 
number of cool PEAGB stars which, as discussed in Sec. 
3.1, can make a significant contribution to the flux of young 
populations. 



3.3 Lick spectral absorption feature indices 

In Table |A"T1 we present all resulting Lick/IDS spectral ab- 
sorption feature indices for Model A' for the seven metallic- 
ities we have considered. In Fig. HOI we show the correspond- 
ing evolutionary curves. We see that all indices except for 
H,3 increase with increasing age and metallicity, with greater 
variation at early ages. For comparison, in Fig. 1101 we also 
give Lick/IDS indices for Model B'. This comparison shows 
that all resulting Lick/IDS indices with binary interactions 
are bluer than those without binary interactions and it ap- 
pears that these changes in Lick indices due to binaries, 
are very small compared to the typical changes with age 
and metallicity. The reason that the inclusion of binary in- 
teractions makes all resulting Lick/IDS spectral absorption 
feature indices bluer is that all of the indices are within the 
wavelength region of 4000 to 6400 A, and in this region the 
ISED for Model A' is bluer than that for Model B' (see Fig. 



4 INFLUENCE OF BINARY INTERACTIONS, 
INPUT PARAMETERS/DISTRIBUTIONS 

We have performed six additional sets of Monte Carlo sim- 
ulations in order to understand the effects of varying the 
model input parameters and initial distributions on the re- 
sulting U-B and B-V colours. Here we focus only on solar 
metallicity BSPs. The parameters we consider are the CE 
ejection efficiency «ce and the Reimers' wind mass-loss co- 
efficient rj. The input distributions we vary are the eccen- 
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Figure 10. The effects of binary interactions on absorption-line indices in the Lick/IDS system for BSPs of various metallicity. The 
symbols linked by a line are for BSPs with binary interactions (Model A'), those without a line are for BSPs without binary interactions 
(Model B'). Different symbol denotes different metallicity, from top to bottom, the metallicity Z is 0.03, 0..004, 0.001, 0.0003 and 0.0001, 
respectively. 
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tricity distribution and the initial mass distribution of the 
secondaries. 

Each model differs from the other by changing one of 
the parameters or distributions used for the initial condi- 
tions of the BSP. The characteristics of each model are sum- 
marized in Table [5] where: the first column is the name 
of each model; the second column denotes the condition 
that binary interactions (Bis) are taken into account; the 
third column gives the eccentricity distribution, where 'E' 
means all binaries are formed in eccentric orbits and the 
eccentricity satisfies a uniform distribution (see equation [fjj 
and 'C denotes initially circular orbits (see equation^; the 
fourth and fifth columns give the values of oce and rj; the 
sixth indicates the initial mass distribution of the secon- 
daries, where 'CC represents the case that secondary-mass 



M2 is correlated with primary-mass Mi and the mass-ratio 
distribution n(q) satisfies a constant form (see equation 
'CQ' represents the correlated and thermal distribution case 
(see equation [3J, and 'UN' is for the uncorrelated case; and 
the last column gives the metallicity Z (solar in all cases) . 

Our standard is Model A which is Model A' restricted 
to solar metallicity. Similarly Model B is Model B' re- 
stricted to only solar metallicity. For Model C binary in- 
teractions are taken into account but all binary orbits are 
initially circular. Model D is the same as Model A except 
that qce = 3.0. Model El does not include any stellar 
wind (77 = 0.0) while in Model E2 77 = 0.5. To investi- 
gate the changes produced by altering the initial secondary- 
mass distribution Model Fl chooses the secondary-mass 
independently from the same IMF as used for the primary 
while in Model F2 the secondary- mass is correlated with 
primary-mass and the mass-ratio satisfies a thermal distri- 
bution. 

To discuss the effects of binary interactions and the in- 
put parameters/distributions on the appearance of BSPs we 
use the difference of the integrated colours (A(U — B)j_a 
and A(B — V)j-a) obtained by subtracting the colour for 
Model A from that of Model j {j =B-F). All results for 
Model A have been given in Section 3 and the differences 
are plotted in Figs II II and 1121 



4.1 Binary interaction (Model B) 

Binary stars play a very important role in EPS studies and 
as we pointed out earlier the majority of current EPS stud- 
ies have only included single star evolution. In Section 3 we 
have already discussed the effect of the modelling of binary 
interactions on the integrated colours, ISED and Lick/IDS 
absorption line indices of a population for a range of metal- 
licities. From the comparison of the integrated U — B and 
B — V colours for Models A and B in Figs ^3 an d Q2I we 
see that the colours for Model A are less than those for 
Model B in almost all instances (the exception is B — V at 
13 Gyr), i.e. for both colours considered the solar metallic- 
ity BSPs including binary interactions are bluer than those 
without binary interactions. For the full metallicity range 
(Z = 0.0001,0.003,0.001,0.004,0.01 and 0.03) Figs □ and H 
show that binary interactions play the same effect on U — B 
and B — V colours. Also Fig. 1101 shows that binary inter- 
actions make all Lick/IDS absorption line indices bluer for 
BSPs regardless of metallicity. The reason that binary inter- 
actions make the integrated U — B and B — V colours and 
the Lick/IDS absorption indices bluer has been discussed in 
Section 3 and we do not take it further here. 



4.2 The distribution of orbital eccentricity 
(Model C) 

Observations of binary stars show that the two compo- 
nents of the binary often move in eccentric orbits and 
that there are some systems for which the eccentricity ap- 
proaches unity. Thus it would be wrong to neglect eccen- 
tricity when mode lling binary evolution and accordingly 
iHurlev et alJ i2002f) included eccentric orbits in their BSE 
algorithm. Statistical studies for the distribution of orbital 
eccentricity of binaries have found that it appears to be 
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Figure 13. The total flux and the the fluxes of all MS stars with 
mass Af >0.7Mq, CHeB and GB stars for Models A (full line) 
and C (dotted line) at an age of 1 Gyr. Each flux is expressed in 
units of the total flux at 2.2 /an (F\ u.2u.m) for the corresponding 
model. 
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Figure 14. Similar to Fig. I13l but for an age of 13 Gyr. 



roughly uniform in the range from zero to unity (see equa- 
tion E| j a t least for those binaries whose components are 
sufficiently well separated that they can have had little 
chance yet to interact with each other, e ven at periastron 
(Eggl eton. Han fc Kiseleva-Eggletonll2004h . So in our stan- 
dard model (Model A) we assumed that all binaries are 
formed with an eccentricity drawn from an uniform eccen- 
tricity distribution. However, to give an idea of how sensitive 



the integrated U — B and B — V colours are to changes in 
this distribution in Model C we consider the rather drastic 
case that all binaries are initially circular initially. 

In Figs II H and 1121 we present the differences in the inte- 
grated U — B and B — V colours between Model C and Model 
A, i.e. A(U-B)( C -A) and A(B-V) ( c-A)- Fig.HTlshows that 
the integrated U — B colour for Model C is greater than that 
for Model A at most ages (except at r = 1,9 and 15 Gyr), 
which means that if all binaries are in circular orbits the 
resulting U — B colour is redder, and thus, the populations 
appear older. In Fig. 1121 a similar effect on the integrated 
B — V colour appears except at r = 1, 9, 13 — 15 Gyr. 

Why does Model C make the integrated U — B and 
B — V colours redder? In the BSE algorithm tidal inter- 
actions are modelled (circularisation and synchronisation) 
which means that an eccentric orbit may be circularised, if 
the stars are close enough for tides to become strong, and 
the orbital separation decreased. Now consider a binary in 
Model C which has a circular orbit and a separation a such 
that mass-transfer is just avoided. The same binary in Model 
A with an initial eccentricity will experience tidal circulari- 
sation (most likely when the primary star reaches the GB) 
which will decrease a and lead to mass-transfer. Also, closer 
systems that would initiate mass-transfer in Model C but 
not merge now are likely to lead to coalescence in Model 
A. In summary, the variation of the eccentricity distribution 
not only can influence the evolutionary path of close binaries 
but also that of relatively wide ones. 

For binary populations, the variation in evolutionary 
path of binaries will cause variations in the number of stars 
and the distribution of stars in the CMD for different evo- 
lutionary phases. These two effects determine the discrep- 
ancies in the integrated colours and ISEDs. Because it is 
extremely difficult to describe the variation of the distribu- 
tion of stars in the CMD we will focus only on the vari- 
ation of the numbers of MS stars with mass M>O.7M0, 
CHeB stars and GB stars. The reason is that the ISED in 
the UBV passbands is mainly contributed by these three 
evolutionary phases for both young (although GB stars do 
not contribute so strongly) and old BSPs (see Figs [3] and 
0J. For r = lGyr BSPs we find that: (i) on the MS the 
number of primaries with mass M > 0.7Mq for Model C is 
less than that for Model A (i.e., AiV Pi MS,(C-A) < 0) while 
that of secondaries with mass M > 0.7Mm for Model C is 
greater (i.e., AJV Si ms (c-A) > 0) and the decrease in the 
number of primaries is less than the increase in the num- 
ber of secondaries (i.e., AJV PiMSi ( C -A) + A -^s,MS,(c-A) > 0); 
(ii) on the GB both the numbers of primaries and secon- 
daries almost do not vary; and (iii) during CHeB the num- 
ber of primaries for Model C is less than that for Model 
A (AAf pi cHeB,(C-A) < 0) while the number of secondaries 
is greater ( A AT S] CHeB, (c- A) > 0) and AAf piCHc B,(C-A) ~ 
AA?s,chcb,(c-a)- For r = 13 Gyr BSPs on the MS and GB 
both the numbers of primaries and secondaries for Model C 
are greater than those for Model A, while the opposite holds 
during CHeB. 

In Figs 1131 and 1141 we give the total flux and the fluxes 
of all MS stars with mass M > O.7M , CHeB stars and GB 
stars over a wavelength range, 4000 < A < 10000A, for Mod- 
els A and C at ages of 1 Gyr and 13 Gyr, respectively. Each 
flux curve is expressed in units of the total flux at 2.2 /im 
(i.e., F\ : 2.2nm) for the corresponding model. Fig. 1131 shows 
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Figure 15. The total flux and the the fluxes of all MS stars with 
mass Af >0.7Mq, CHeB and GB stars for Models A (full line) 
and D (dotted line) at an age of 1 Gyr. Each flux is expressed in 
units of the total flux at 2.2 /an (i*A,2.2/jm) f°r the corresponding 
model. 



that the total flux and the fluxes eminating from the three 
evolutionary phases are redder for Model C than for Model 
A at an age of r = 1 Gyr. Comparing the discrepancies in the 
fluxes between Models C and Model A we can conclude that 
the redder total flux (the top panel of Fig. 1131 and the redder 
colours (Figs lllj and ll21 for Model C are mainly introduced 
by the difference in MS stars with mass M > 0.7Mq at an 
age of r = lGyr. For the r = 13 Gyr BSP.~Fig. [HI shows 
that the total flux and the flux of all CHeB stars for Model 
C are redder than for Model A, and the discrepancies in the 
fluxes of all MS stars with mass M > O.7M0 and GB stars is 
insignificant. Thus it is the CHeB stars that are influencing 
the redder total flux (top panel of Fig. 1141 and the redder 
colours (Figs [TT] and [TU for Model C at 13 Gyr. 

4.3 The common-envelope ejection efficiency oce 
(Model D) 

Common-envelope evolution is one of the most important 
and complex but also one of the least understood phases 
of binary evolution. One area of the uncertainties is related 
to the criterion for the ejection of the CE, which crucially 
dete rmines the orbit al period distribut ion of post-CE bina - 
ries <Han et al.ll2002t . In the BSE code lHurlev et al.l (l2002t 
adopt a widely used and relatively simple criterion where 
the CE is ejected when the change in orbital energy, multi- 
plied by the CE ejection efficiency parameter, lvce, exceeds 
the binding ene rgy of the envelop e (see equation|7|l. The ap- 
proach used bv lHan et alJ <|2002T) differs in that the thermal 
energy of the envelope is also considered in their criterion. 
The value of cxce is an uncertain but very crucial factor as it 
determines the evolution path of post-CE binaries and there- 
fore affects the birth rates and numbers of the various types 
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Figure 16. Similar to Fig. 1151 but for an age of 13 Gyr. 



of binary in binary population synthesis (B P S). F or exam- 
ple, using the BSE algorithm iHurlev et alJ <l2002l) showed 
that a particular post-CE binary that went on to form a 
cataclysmic variable when oce = 1 was used would emerge 
from CE evolution with too wide an orbit for subsequent 
interaction if «ce = 3 was used instead. Knowledge of the 
typical values of «ce is therefore crucial in understanding 
the evolution of populations of binary systems. Comparing 
observational birth rates and numbers of the various binary 
types with those obtained by the method of BPS, we can 
set constrains on oce theoretically, but in fact this method 
is difficult because selection effects involved in the obser- 
vations weaken the constraints. In this work we vary oce 
in a reasonable range to investigate its effect: in Model A 
qce = 1, while in Model D it is set to 3. 

Figs llll and ll2| give the difference in the integrated U— B 
and B — V colours between Models D A, i.e., A(U — B)(d-a) 
and A(B — V)m_A)- Fig. Illl shows that the integrated U — 
B colour for Model D is less than that for Model A in all 
instances and Fig. ll2l shows that the integrated B — V colour 
for Model C is less than that for Model A except at r = 2 and 
8 Gyr. That is to say, by increasing oce the resulting U — B 
and B — V colours have become bluer and the population 
appears younger. The effect of Model D on the integrated 
U — B and B — V colours is contrary to that of Model C. 

Because cecE is introduced in CE evolution its varia- 
tion only affects the evolution of binaries experiencing a 
CE phase: it affects the CE evolutionary process and the 
evolution path of post-CE descendants. If qce is high less 
energy is required to drive off the CE and thus it is eas- 
ier to eject the CE. This means that CE ejection occurs 
earlier in the interaction with the result that most of the 
post-CE binaries will systematically have relatively longer 
orbital periods. However, if qqe is much lower than unity, 
then short-period CE descendants along with mergers will 
be much more common. In Model D with acE = 3 the or- 
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bital periods of post-CE binaries systematically increases so 
that close pre-CE binary systems that would have coalesced 
during CE in Model A now survive and go on to interact. 
On the other hand, post-CE binary systems that would lead 
to mass-transfer in Model A now go through their entire 
evolution independently. This variation in the evolutionary 
path of post-CE binaries causes variations in the number of 
stars and their distribution in the CMD for different evolu- 
tionary phases and this leads to the observed discrepancies 
in the integrated colours and ISEDs. We find that overall 
the number of binaries experiencing RLOF does not change 
when qce increases but that the number experiencing more 
than one phase of RLOF increases significantly. 

In Figs 1151 and 1161 we give the total flux, and the fluxes 
of all MS stars with mass M>O.7M , CHeB stars and 
GB stars over the wavelength range 4000 < A < 10000A 
for Models A and D at ages of 1 Gyr and 13 Gyr, respec- 
tively. Once again each flux curve is expressed in units of 
the total flux at 2.2 /im (i.e., F\ t 2.2pm) for the correspond- 
ing model. Fig. 1151 shows that the flux of all MS stars with 
mass M > 0.7M o for Model D is bluer than that for Model 
A, and there are no significant deviations in the GB and 
CHeB fluxes. So the bluer ISEDs and integrated (U - B) 
and (B - V) colours (see Figs [TT] and for Model D at 
1 Gyr are mainly the result of the differences in MS stars 
with mass M > O.7M . For the r = 13 Gyr BSPs, Fig. fTTH 
shows that the flux of all MS stars with mass M > O.7M0 
for Model D is also bluer than that of Model A and the 
CHeB flux is redder. The net effect is the bluer ISEDs and 
integrated (U — B) and (B — V) colours. 



4.4 Reimers' wind mass-loss efficiency (Model 
El-2) 

Mass-loss is one of the most important processes in stel- 
lar astrophysics. Its value often determines the evolutionary 
path of a star and therefore the final stages of a binary sys - 
tem. The BSE/SSE packages of iHurlev et~al"l feOOOt I2OO2I) 
included a prescription for mass lo s s whic h is not included in 
the detailed models of iPols et alJ (ll99ST) . This prescription 
is drawn from a range of current mass-loss theories available 
in the literature and can easily be altered or added to (for 
details refer to Hurley et al. 2000). 

We on ly discuss the in fluence of Reimers' empirical wind 
mass- loss jReimers!ll975T) on our results because it is the 
most relevant among these mass-loss mechanisms for pop- 
ulations of age t > 1.0 Gyr. Reimers' wind mass-loss has 
been widely used in many EPS studies owing to its simple 
parameterized form, and in the SSE/BSE packages it is ap- 
plied to the stellar envelope for intermediate- and low-mass 
stars on the GB and beyond. The value of the Reimers' 
mass- loss efficiency fact or 77 has been quoted previously as 
varying from 0.25 to 2-3 <DuDreJl986tlKudritzki fc Reimers! 
ll97Sl: lR cnzm ] l!98lft . The possibility of a metallicity depen- 
dence it is quite controversial: it has been argued that there 
is a positive metallicity dependen ce supported by observa- 
tions and hydrodynamical models iYi. Demaraue fc Oemlerl 
Il997l) . while others have stated that a metallicity depen- 
dence should not be included as there is no strong evidence 
for it llben fc Renzinil Il983t ICarraro et al.lll99a) . In this 
study we do not include a metallicity dependence in the 



Reimers' mass- loss, which is to say that we adopt a fixed 
r) across all metallicities in any particular model. In order 
to investigate the effects of wind-loss on the model colours, 
three reasonable choices for the Reimers' coefficient are used 
in equation ©: T)= 0.0 (Model El), 0.3 (Model A), and 0.5 
(Model E2). 

In Figs 1111 and 1121 the differences in the integrated 
U — B and B — V colours between Model El and Model 
A (i.e., A(U - B) (E1 _ A) and A(B - V)( E i-A)) and tnose 
between Model E2 and Model A (i.e., A(U - B) (E2 _ A ) and 
A(B — V)(E2-A)) are given. By comparing the integrated 
U — B and B — V colours for Models El and A we see 
that these two colours for Model El fluctuate around those 
for Model A at early age (r < 8 Gyr), at intermediate and 
late ages the values of these two colours for Model El are 
larger/redder than those for Model A, and the discrepancies 
at intermediate and late ages are greater than those at early 
age. Comparison of the integrated U — B and B — V colours 
between Models E2 and A shows that there is no systematic 
discrepancy in the two colours. 

In summary, if neglecting the Reimers' mass-loss alto- 
gether (Model El) the integrated U — B and B — V colours 
are redder and the populations look older at intermediate 
and late ages, while at early age there is no systematic dis- 
crepancy (r < 8 Gyr). The variation of 77 from 0.3 to 0.5 
(Model E2) does not give rise to systematic differences in 
the integrated U — B and B — V colours. 

The reason that the variation of r\ can influence the in- 
tegrated colours, is that the variation in mass-loss owing to 
a stellar wind can influence the spin orbital angular momen- 
tum of a star and also the total orbital angular momentum. 
This is in addition to the obvious effect it has on the stellar 
mass as well as possibly the companion mass via accretion. 
So variation in 77 can alter the appearance of a star and even 
the evolution path of a binary. This in turn will lead to a 
variation in the distribution of stars in the CMD and the 
integrated colours, especially in the extreme case of 77 = 
(Model El). 

4.5 The initial mass distribution of the 
secondaries (Model Fl-2) 

The distribution of mass-ratio is less well-known than either 
of the distributions over period or stellar mass. This is be- 
cause substantially more orbit data is required to determine 
a mass ratio than for the total mass or period of a binary. 
Therefore, it is quite controversial. In Model Fl we assume 
that the primary- and secondary-ma sses are uncorrelated 
(both are chosen from equation 0- iDuauennov fc Mavorl 
( 1991) found this to be an adequate approximation for their 
sample of binaries whose primaries were all F/G dwarfs like 
the Sun but it cannot be an adequat e approximation for 
massive stars in short-period (P < 25 d. lLucv fc Riccc^ l979) 
and moderately-shor t-period (P < 3000 d, Maze h et alJ 
Il992t lTokovininiri992l) binaries. In Models F2 and A we as- 
sume that the two component masses are correlated: Model 
A takes a constant mass-ratio distribution (see equation |2J 
while Model F2 takes a thermal form (see equation [3J . 

From the discrepancy in the integrated (U — B) and 
(B-V) colours between Model Fl and Model A (i.e., A(U- 
B)( F i_A) and A(B — V)(fi_a)) and those between Model F2 
and Model A (i.e., A(U - B) (F2 -A) and A(B - V) (F2 -A)) 



14 F. Zhang et al. 



EFT 

uncorrelatecL 

-- n(q) = 1 




-0.5 





log M 



0.5 



Figure 17. The initial mass distributions of primaries (thick) and 
secondaries (thin) for three assumptions: uncorrelated compo- 
nent masses (Model Fl), unique mass-ratio distribution n(q) = 1 
(Model A) and thermal mass-ratio distribution n(q) = 2q (Model 
F2), and the distribution of EFT (1989). Primary- and secondary- 
mass is in the range 0.1 < M < IOOMq. 



in Figs 1111 and 1121 we see that the variation in the initial 
secondary-mass distribution leads to fluctuations in the in- 
tegrated (U — B) and (B — V) colours and that the fluctuation 
at late age is greater than that at early age. 

In Fig. 1171 we present the initial mass distributions of 
primaries and secondaries for Models A, Fl and F2. From it 
we see that massive stars are more likely to have low-mass 
companions for the distribution of uncorrelated component 
masses (Model Fl) than for those correlated cases (Models A 
and F2). Furthermore, the possibility of massive stars paired 
with low-mass stars is greater for a constant (Model A) than 
for a thermal initial mass- ratio distribution (Model F2). In 
Fig. ^3 we a l so present the distribution of EFT (1989, refer 
to equation^ . The difference in the initial distribution leads 
to the discrepancy observed in the integrated (U — B) and 
(B — V) colours. 



5 SUMMARY AND CONCLUSIONS 

We have simulated realistic stellar populations composed of 
100% binaries by producing 1 x 10 6 binary systems using 
a Monte Carlo technique. Using the EPS method we com- 
puted the integrated colours, ISEDs and Lick/IDS absorp- 
tion feature indices for an extensive set of instantaneous 
burst BSPs with and without binary interactions over a 
large range of age and metallicity: 1 Gyr < r < 15 Gyr 
and -2.3 < [Fe/H] < 0.2. In our EPS models we adopted 
the rapid SSE and BSE algorithms for the single and bi- 
nary evolution tracks, the empirical and semi-empirical cal- 
ibrated BaSeL-2.0 model for the library of stellar spectra 
and empirical fitting functions for the Lick/IDS spectral ab- 
sorption feature indices. By comparing the results for pop- 
ulations with and without binary interactions we show that 
the inclusion of binary interactions makes each quantity we 
have considered bluer and the appearance of the population 
younger or the metallicity lower. 



Also we have investigated the effect of input parame- 
ters (the CE ejection efficiency acE and the Reimers' wind 
mass-loss coefficient rf) and the input distributions (eccen- 
tricity and the initial mass of the secondaries) on the in- 
tegrated U — B and B — V colours of a solar-metallicity 
BSP. The results reveal that the variations in these pa- 
rameters/distributions can significantly affect the results. 
Comparing the discrepancies in the integrated colours for 
all models, we find that the differences between the models 
with and without binary interactions are greater than those 
caused by the variations in the choice of input parameters 
and distributions. Based on the above results an important 
conclusion can be drawn that it is very necessary to include 
binary interactions in EPS studies. 

The integrated colours and Lick/IDS absorption-line in- 
dices for Models B-F, and the ISEDs for all models are not 
given for the sake of the length of the paper but are available 
on request. 
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APPENDIX A: THE ABSORPTION INDICES FOR MODEL A 



Ago 


1 


2 


3 


4 


5 ( 


> 7 8 


9 10 


11 


12 


13 


14 


15 


(Gyr) 



























z = 0.0001 



CNi 


-0 


.228 


-0 


.201 


-0 


.178 


-0 


.162 


-0 


-163 


-0 


.152 


-0 


.152 


-0.149 


-0 


.146 


-0 


.154 


-0 


.147 


-0 


.148 


-0 


.149 


-0.146 


-0 


.147 


CN 2 


-0 


.129 


-0 


.116 


-0 


.102 


-0 


095 


-0 


.100 


-0 


.093 


-0 


.094 


-0.091 


-0 


089 


-0 


096 


-0 


.091 


-0 


.092 


-0 


.092 


-0.089 


-0 


.091 


Ca4227 


-0 


.346 


-0 


.215 


-0 


133 


-0 


052 


-0 


.024 


o 


.014 





.022 


0.034 





.071 





068 


o 


.094 


o 


.116 


o 


.111 


0.133 





124 


G4300 


-3 


.215 


-2 


.326 


-1 


.604 


-1 


.097 


-1 


.085 


-0 


.758 


-0 


.752 


-0.648 


-0 


.572 


-0 


.756 


-0 


.558 


-0 


.545 


-0 


.575 


-0.468 


-0 


.478 


Fc4383 


-2 


.066 


-1 


.555 


-1 


082 


-0 


821 


-0 


.819 


-0 


.611 


-0 


.522 


-0.421 


-0 


.399 


-0 


.506 


-0 


.362 


-0 


.353 


-0 


.340 


-0.280 


-0 


.233 


Ca4455 


-0 


.036 





026 





082 





109 


o 


.092 


o 


.117 





.122 


0.137 





.145 





120 


o 


.143 


o 


.144 


() 


.144 


.151 





.147 


Fe4531 


-0 


.290 


-0 


.014 





195 





304 


o 


.241 


o 


.334 





.345 


0.386 





.407 





.308 


o 


.391 


o 


.383 


o 


.388 


0.395 





366 


Fc4668 


1 


.317 





878 





667 





488 


o 


.379 


o 


.250 





.262 


0.257 





.265 





.245 


o 


.222 


o 


.193 


o 


.230 


0.226 





217 


H« 
ll p 


5 


.888 


4 


963 


4 


320 


3 


896 


3 


.948 


3 


.706 


3 


.682 


3.583 


3 


485 


3 


639 


3 


.426 


3 


.393 


3 


.370 


3.259 


3 


.262 


Fe5015 


o 


.606 





812 





.980 


1 


028 


o 


.905 


o 


.975 





.980 


1.015 


1 


026 





.904 


o 


.984 


o 


.971 


o 


.977 


0.971 


{1 


930 


Mgi 





.011 





005 





003 





002 





.001 





.000 





.001 


0.001 





003 





004 





.004 





.005 





.007 


0.008 





009 


Mg 2 


o 


.016 





022 





027 





028 


o 


.029 


o 


.031 





.033 


0.034 





036 





036 


o 


.037 


o 


.039 


o 


.039 


0.040 





041 


iV -"&b 


o 


.373 


(1 


495 





525 


{1 


559 


o 


.731 


o 


.769 





.790 


0.856 


(1 


.820 





.871 


o 


.851 


o 


.871 


o 


.875 


0.887 





907 


Fc5270 


-0 


.807 


-0 


.494 


-0 


270 


-0 


.158 


-0 


.253 


-0 


.159 


-0 


.149 


-0.098 


-0 


075 


-0 


171 


-0 


.073 


-0 


.064 


-0 


.055 


-0.035 


-0 


059 


Fc5335 





.366 





428 





492 





513 


o 


.461 


o 


.487 





.505 


0.530 





.545 





.507 


o 


.544 


o 


.553 


o 


.559 


0.570 





561 


Fc5406 


o 


.107 





148 





195 





205 


o 


.153 


o 


.173 





.184 


0.204 





215 





176 


o 


.209 


o 


.212 


o 


.219 


0.224 





214 


Fc5709 


o 


.094 





.104 


(1 


125 


(1 


132 


o 


.105 


o 


. 1 1 1 





.110 


0.123 





122 





102 


o 


.119 


o 


.114 


o 


.120 


0.122 





115 


Fc5782 


-0 


.217 


-0 


.187 


-0 


158 


-0 


148 


-0 


.180 


-0 


.169 


-0 


.165 


-0.153 


-0 


153 


-0 


180 


-0 


.159 


-0 


.163 


-0 


.160 


-0.160 


-0 


169 


Na D 


1 


.149 


1 


050 


1 


095 


1 


002 


o 


.853 


o 


.886 





.911 


0.934 





.951 





902 


o 


.932 


o 


.999 


o 


.997 


1.001 





997 


TiOi 


o 


.014 





013 





013 





012 


o 


.012 


o 


.011 





.012 


0.011 





Oil 





.011 


o 


.010 


o 


.011 


o 


.010 


0.010 





010 


Ti02 


o 


.000 





000 





000 


-0 


001 


-0 


.003 


-0 


.003 


-0 


.002 


-0.004 


-0 


003 


-0 


004 


-0 


.004 


-0 


.004 


-0 


.004 


-0.005 


-0 


005 




























Z = 


0003 




























CNi 


-0 


.205 


-0 


.181 


-0 


.156 


-0 


.137 


-0 


.125 


-0 


.111 


-0 


.108 


-0.105 


-0 


097 


-0 


099 


-0 


.096 


-0 


.101 


-0 


.098 


-0.106 


-0 


.107 


CN 2 


-0 


.117 


-0 


.103 


-0 


.087 


-0 


.075 


-0 


.070 


-0 


.062 


-0 


.060 


-0.060 


-0 


053 


-0 


052 


-0 


.052 


-0 


.055 


-0 


.051 


-0.058 


-0 


058 


Ca4227 


-0 


.191 


-0 


.062 





033 





094 


o 


.142 


o 


.194 


o 


.199 


0.203 


(1 


265 


(1 


.266 


o 
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o 
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-2 


.248 


-1 


.395 


-0 
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-0 


.011 


o 


.318 


o 


.775 





.902 


0.998 


1 


.235 


1 


.139 


1 


.236 


1 


.157 


1 


.205 
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.975 


Fc4383 


-1 


.431 


-1 


.081 


-0 


672 


-0 


287 


-0 


.152 


o 
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.306 


0.319 





.528 
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o 


.640 
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o 


.689 


0.572 





.571 


Ca4455 


-0 


.058 





065 





.148 





217 
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.269 





.280 


0.272 





310 





.317 


o 


.323 


o 
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o 


.345 


0.309 





.311 


Fc4531 


o 
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.468 
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861 


o 
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o 
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o 
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.240 
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.260 
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.259 


-0 
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-0 


.286 


-0 


.224 


-0.270 


-0 
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5 


.495 


4 


.776 


4 


152 


3 


.702 


3 


.479 


3 


.140 


3 


.013 


2.945 


2 


710 


2 


740 


2 


.664 


2 


.727 


2 


.657 


2.786 


2 


770 


Fe5015 





.735 


1 


146 


1 


383 


1 


542 


1 


.476 


1 


.571 


1 


.599 


1.499 


1 


626 


1 


676 


1 


.651 


1 


.636 
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.757 


1.609 


1 
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o 
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o 
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.005 


o 


.006 


o 


.007 
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.009 


0.009 


() 
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o 


.029 





037 





042 





.047 
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.048 
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050 
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o 
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o 


.056 
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o 


.573 


(1 


635 





725 





904 


o 


.886 


o 


.950 





.958 


1.055 


1 


022 


1 


019 


1 


.069 


1 


.086 


1 


.058 


1.045 


1 


058 


Fe5270 


-0 


.330 
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265 





.420 


o 


.411 


o 


.499 





.528 


0.477 





576 





612 


o 


.608 


o 


.606 


o 


.693 


0.606 





605 


Fc5335 





.273 





422 





524 





597 


o 


.554 


o 


.603 





.620 


0.585 





.640 





661 


o 


.667 


o 


.678 


o 


.734 


0.685 





683 


Fc5406 


o 


.087 





173 





240 





.281 


o 


.234 


o 


.263 





.273 


0.236 





280 





303 


o 


.301 


o 


.308 


o 


.356 


0.314 





309 


Fe5709 


o 


.127 





173 





216 





243 


o 


.243 


o 


.252 





.253 


0.238 





260 





.277 


o 


.270 


o 


.261 


o 


.288 


0.271 


{1 


.269 


Fc5782 


-0 


.075 


-0 


.017 


(1 


017 





046 


o 


.026 


o 


.040 





.043 


0.021 


(1 


.044 





058 


o 


.052 


o 


.048 
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.075 


0.051 





045 
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1 


.052 


1 


067 


1 


086 


1 


095 


o 


.874 


o 


.911 





.949 
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.896 





962 


o 


.917 


1 


.018 


1 


.118 


1.007 





.989 


TiOi 


o 


.018 





020 





019 





018 


o 


.014 


o 


.015 





.015 


0.013 





.013 





013 


o 


.013 


o 


.014 


o 


.014 


0.013 


o 


012 
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o 


.009 





014 





013 





013 


o 
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o 


.005 





.006 


0.001 





003 





002 


o 


.002 


o 


.004 


o 


.006 


0.002 





001 
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001 




























CNi 


-0 


.192 


-0 


.172 


-0 


136 


-0 


.118 


-0 


.101 


-0 


.088 


-0 


.081 


-0.075 


-0 


076 


-0 


076 


-0 


.071 


-0 


.082 


-0 


.086 


-0.086 


-0 


098 


CN 2 


-0 


.110 


-0 


.097 


-0 


.075 


-0 


062 


-0 


.051 


-0 


.043 


-0 


.039 


-0.035 


-0 


037 


-0 


038 


-0 


.032 


-0 


.041 


-0 


.045 


-0.045 


-0 


051 


Ca4227 


-0 


.039 





123 





251 





388 


o 


.403 


o 


.395 





.446 


0.491 





.443 


(1 


.442 


o 


.475 


o 


.445 


o 


.377 


0.457 


{1 


392 


G4300 


_1 


.488 


-0 


.791 





417 


1 


122 


1 


.654 


2 


.053 


2 


.294 


2.488 


2 


.468 


2 


.454 


2 


.599 


2 


.283 


2 


.161 


2.208 


1 


.759 


Fc4383 


-0 


.955 


-0 


.660 


-0 


003 





494 


o 


.828 


1 


.136 


1 


.209 


1.365 


1 


341 


1 


332 


1 


.422 


1 


.416 


1 


.341 


1.380 


1 


341 


Ca4455 


o 


.017 





169 


o 


315 





422 


o 


.491 


o 


.511 





.533 


0.569 





545 





.540 


o 


.581 


o 


.526 


o 


.510 


0.529 





503 


Fc4531 


o 


.522 


(1 


.880 


1 


.234 


1 


.448 


1 


.563 


1 


.624 


1 


.657 


1.731 


1 


673 


1 


666 


1 


.762 


1 


.630 


1 


.605 


1.643 


1 


.606 


Fc4668 


-0 


.356 


-0 


.148 


-0 


171 





084 


o 


. 143 


-0 


.015 





.017 


0.093 


-0 


019 


-0 


058 


o 


.005 


-0 


.082 


-0 


.156 


-0.158 


-0 


081 


i_i_p 


5 


.269 


4 


684 


3 


893 


3 


.415 


3 


.080 


2 


.797 


2 


.629 


2.488 


2 


496 


2 


.474 


2 


.369 


2 


.519 


2 


.639 


2.585 


2 


834 


Fe5015 


1 


.135 


1 


810 


2 


055 


2 


300 


2 


.493 


2 


.406 


2 


.408 


2.506 


2 


395 


2 


365 


2 


.518 


2 


.315 


2 


.261 


2.312 


2 


329 


Mgx 


o 


.015 





014 





013 





016 


o 


.015 


o 


.015 





.016 


(J. 018 





017 





.017 


o 


.019 


o 


.018 


o 


.017 


0.020 


(1 


.021 


Mg 2 





.046 





056 





062 





.075 





.076 





.074 





.076 


0.083 





079 





078 





.084 





.078 





.082 


0.080 





083 


Mg b 





.798 


1 


106 


1 


176 


1 


388 


1 


.574 


1 


.391 


1 


.442 


1.497 


1 


.447 


1 


.417 


1 


.403 


1 


.447 


1 


.616 


1.448 


1 


464 


Fe5270 





.266 





.578 





.821 





991 


1 


.110 


1 


.142 


1 


.167 


1.232 


1 


182 


1 


.178 


1 


.268 


1 


.147 


1 


.125 


1.174 


1 


153 


Fe5335 





.410 





611 





726 





917 





.934 





.935 





.943 


1.024 





937 





935 


1 


.003 





.918 





.891 


0.926 





946 


Fe5406 





.199 





318 





382 





.497 





.502 





.511 





.510 


0.580 





508 





.504 





.563 





.485 





.466 


0.492 





516 


Fe5709 





.230 





274 





338 





337 





.381 





.402 





.411 


0.423 





.415 





.420 





.447 





.408 





.406 


0.416 





.412 


Fe5782 





.099 





133 





178 





212 





.230 





.238 





.241 


0.261 





239 





238 





.269 





.229 





.224 


0.234 





245 


Na D 


1 


.169 


1 


.147 


1 


086 


1 


283 


1 


.224 


1 


.159 


1 


.156 


1.268 


1 


.190 


1 


142 


1 


.266 


1 


.167 


1 


.156 


1.255 


1 


286 


TiOi 





.026 





035 





029 





037 





.031 





.025 





.023 


0.025 





.021 





.019 





.020 





.019 





.017 


0.017 





021 


Ti0 2 





.025 





.042 





031 





.048 





.037 





.025 





.021 


0.026 





017 





.013 





.018 





.014 





.010 


0.011 





.019 






























Z = 


0.004 




























CNi 


-0 


.185 


-0 


.136 


-0 


.105 


-0 


079 


-0 


.064 


-0 


.050 


-0.047 


-0.051 


-0 


038 


-0 


033 


-0 


.031 


-0 


.044 


-0 


.037 


-0.035 


-0 


035 


CN 2 


-0 


.111 


-0 


.077 


-0 


055 


-0 


036 


-0 


.025 


-0 


.014 


-0.011 


-0.014 


-0 


.005 


-0 


.001 


-0 


.001 


-0 


.010 


-0 


.005 


-0.004 


-0 


005 


Ca4227 





.159 





339 





427 





518 





.604 





.652 


0.696 


0.690 





.710 





.792 





.789 





.775 





.778 


0.856 





803 


G4300 


-0 


.955 





661 


1 


690 


2 


564 


3 


.081 


3 


.502 


3.658 


3.581 


4 


.012 


4 


258 


4 


.316 


3 


.978 


4 


.207 


4.366 


4 


359 


Fe4383 


-0 


.406 





375 


1 


101 


1 


677 


2 


.017 


2 


.322 


2.429 


2.482 


2 


684 


2 


.873 


2 


.878 


2 


.820 


3 


.020 


3.098 


3 


156 


Ca4455 





.212 





.478 





630 





.754 





.852 





.919 


0.961 


0.958 





983 


1 


037 


1 


.058 


1 


.048 


1 


.059 


1.103 


1 


.073 


Fc4531 


1 


.040 


1 


.553 


1 


.825 


2 


055 


2 


.202 


2 


.327 


2.379 


2.352 


2 


.429 


2 


.524 


2 


.536 


2 


.501 


2 


.558 


2.597 


2 


.570 


Fc4668 


-0 


.096 





409 





800 


1 


030 


1 


.355 


1 


.460 


1.639 


1.703 


1 


506 


1 


592 


1 


.716 


1 


.822 


1 


.701 


1.891 


1 


678 
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Table Al: continued 



Ago 
(Gyr) 




1 




2 




3 




4 




5 




6 




7 




8 




9 




10 




11 




1 2 




13 




14 




15 


H 


5 


.502 


4 


338 


3 


563 


2 


963 


2 


.670 


2 


.398 


2 


339 


2 


384 


2 


147 


2 


037 


2 


.014 


2 


.174 


2 


.018 


1 


965 


1 


945 


Fe5015 


2 


.102 


2 


967 


3 


313 


3 


473 


3 


.815 


3 


.878 


4 


.098 


4 


168 


3 


.874 


3 


977 


4 


.146 


4 


.300 


3 


.996 


4 


243 


3 


952 


Mgl 





.012 





014 





019 





025 





.027 





.032 





.033 





033 





.038 





042 





.042 





.042 





.047 





048 





048 


Mg 2 





.064 





080 





092 





.101 





.112 





.117 





.123 





124 





124 





133 





136 





.140 





146 





147 





.144 


Mg b 


1 


.211 


1 


552 


1 


754 


1 


.844 


2 


.134 


2 


.084 


2 


.224 


2 


306 


2 


127 


2 


232 


2 


.358 


2 


.522 


2 


.530 


2 


628 


2 


471 


Fe5270 





.829 


1 


246 


1 


480 


1 


652 


1 


.789 


1 


.886 


1 


.953 


1 


940 


1 


963 


2 


035 


2 


.060 


2 


.076 


2 


.088 


2 


143 


2 


110 


Fe5335 





.641 





922 


1 


105 


1 


243 


1 


.364 


1 


.441 


1 


.501 


1 


507 


1 


507 


1 


582 


1 


.604 


1 


.634 


1 


.662 


1 


702 


1 


666 


Fe5406 





.296 





492 





617 





723 





.793 





.868 





.888 





883 





918 





969 





.971 





.974 


1 


.020 


1 


027 


1 


013 


Fe5709 





321 





450 





516 





582 





.591 





.643 





.630 





608 





668 





675 





.668 





.645 





.683 





671 





685 


Fo5782 





.214 





299 





353 





406 





.422 





.462 





.461 





.447 





.475 





495 





.486 





.476 





.509 





499 





502 


Na D 


1 


.235 


1 


233 


1 


304 


1 


374 


1 


.501 


1 


.510 


1 


.615 


1 


651 


1 


598 


1 


769 


1 


726 


1 


.802 


1 


.843 


1 


893 


1 


893 


TiOi 





.035 





035 





035 





031 





037 





.033 





.039 





043 





.031 





035 





.038 





.043 





035 





039 





033 


Ti0 2 





.041 





042 





041 





036 





.047 





.041 





.053 





059 





038 





046 





.051 





.061 





.047 





054 





042 






























Z — 


0.01 
































CNx 


-0 


.155 


-0 


.087 


-0 


.056 


-0 


045 


-0 


.036 


-0 


.021 


-0 


.022 


-0 


.012 


-0 


.005 


-0 


006 


-0 


.003 


-0 


002 





001 





006 





004 


CN 2 


-0 


.092 


-0 


.041 


-0 


017 


-0 


007 





.000 


o 


.013 


o 


.014 


o 


022 


o 


028 





027 





028 





030 





032 





035 





034 


Ca4227 





.276 


o 


518 





644 


o 


695 





.720 


o 


.820 


o 


.826 


o 


.875 


o 


919 





942 





.948 





.955 





.968 


1 


.007 





992 


G4300 





.051 


2 


221 


3 


186 


3 


524 


3 


.800 


4 


.221 


4 


.179 


4 


450 


4 


672 


4 


707 


4 


.808 


4 


820 


4 


.913 


5 


121 


5 


(176 


Fc4383 





.379 


1 


805 


2 


575 


2 


912 


3 


.302 


3 


.579 


3 


.647 


3 


839 


4 


027 


4 


082 


4 


.149 


4 


.228 


4 


.313 


4 


400 


4 


422 


Ca4455 


o 


.463 





828 





994 


1 


067 


1 


127 


1 


.213 


1 


.218 


1 


278 


1 


.325 


1 


340 


1 


.352 


1 


.364 


1 


.389 


1 


421 


1 


425 


Fe4531 


1 


.496 


2 


148 


2 


430 


2 


545 


2 


.646 


2 


.784 


2 


.783 


2 


885 


2 


959 


2 


977 


2 


986 


3 


010 


3 


.048 


3 


076 


3 


091 


Fe4668 





.613 


1 


795 


2 


304 


2 


544 


2 


.711 


2 


.971 


2 


.985 


3 


137 


3 


252 


3 


279 


3 


.290 


3 


.291 


3 


.339 


3 


415 


3 


382 


H /3 


5 


.133 


3 


.487 


2 


836 


2 


621 


2 


.475 


2 


.254 


2 


.270 


2 


109 


1 


995 


1 


976 


1 


.930 


1 


.897 


1 


.854 


1 


.774 


1 


767 


Fo5015 


2 


.768 


3 


745 


4 


070 


4 


245 


4 


.342 


4 


.494 


4 


.508 


4 


613 


4 


686 


4 


733 


4 


.750 


4 


.730 


4 


.784 


4 


874 


4 


863 


Mgi 





.014 





026 





037 





041 





.047 





.053 





.054 





060 





064 





065 





.065 





.069 





.070 





071 





074 


Mg 2 





.076 





111 





129 





138 





.147 





.159 





.162 





169 





177 





181 





182 





.185 





189 





.191 





194 


Mg b 


1 


.323 


1 


961 


2 


137 


2 


299 


2 


.388 


2 


.546 


2 


.592 


2 


676 


2 


764 


2 


850 


2 


.896 


2 


.905 


2 


.952 


3 


042 


3 


020 


Fo5270 


1 


262 


1 


815 


2 


053 


2 


160 


2 


.247 


2 


.356 


2 


.372 


2 


449 


2 


507 


2 


534 


2 


.547 


2 


.565 


2 


.595 


2 


622 


2 


639 


Fo5335 





.963 


1 


450 


1 


669 


1 


765 


1 


850 


1 


.956 


1 


.962 


2 


041 


2 


099 


2 


122 


2 


.127 


2 


.150 


2 


.178 


2 


193 


2 


213 


Fo5406 





.510 





840 


1 


009 


1 


079 


1 


146 


1 


.228 


1 


.237 


1 


297 


1 


343 


1 


348 


1 


.352 


1 


.378 


1 


.394 


1 


396 


1 


414 


Fe5709 





.513 





655 





737 





769 





.794 





.825 





.828 





853 





869 





862 





.865 





.881 





.884 





.877 





889 


Fe5782 





.346 





491 





561 





590 





.616 





.649 





.652 





676 





692 





689 





.683 





.695 





.700 





689 





699 


Na D 


1 


.256 


1 


606 


1 


789 


1 


879 


1 


.976 


2 


.086 


2 


.153 


2 


186 


2 


257 


2 


338 


2 


.347 


2 


.355 


2 


.390 


2 


431 


2 


.457 


TiOi 





.027 





031 





030 





031 





.032 





.034 





.034 





034 





034 





036 





.035 





.033 





.034 





036 





036 


TiQ 2 





.024 





035 





036 





038 





.040 





.044 





.044 





045 





047 





050 





.047 





.044 





.047 





049 





049 



Z = 0.02 



CNi 


-0 


121 


-0 


061 


-0 


030 


-0 


011 





.004 





.015 





027 





020 





033 





046 





048 





.042 





056 





030 





063 


CN 2 


-0 


.065 


-0 


.018 





008 





025 





.038 





.048 





.060 





055 





065 





077 





.080 





.074 





.089 





068 





095 


Ca4227 





.411 





690 





834 





917 





.994 


1 


.030 


1 


.126 
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Table Al: continued 



Ago 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


(Gyr) 
































TiOi 


0.024 


0.024 


0.029 


0.033 


0.032 


0.035 


0.034 


0.040 


0.043 


0.037 


0.039 


0.039 
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0.041 
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TiQ 2 
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0.027 


0.038 


0.048 


0.047 


0.053 
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0.065 


0.070 


0.060 


0.064 


0.065 


0.067 


0.068 


0.070 
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